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Abstract 

Collective flow and femtoscopy in ultrarelativistic 3 He-Au collisions are investigated within the 3+1-dimensional (3+1D) 
viscous event-by-event hydrodynamics. We evaluate elliptic and triangular flow coefficients as functions of the transverse 
momentum. We find the typical long-range ridge structures in the two-particle correlations in the relative azimuth and 
pseudorapidity, in the pseudorapidity directions of both Au and 3 He. We also make predictions for the pionic interfero¬ 
metric radii, which decrease with the transverse momentum of the pion pair. All features found hint on collectivity of 
the dynamics of the system formed in 3 He-Au collisions, with hydrodynamics leading to quantitative agreement with 
the up-to-now released data. 

Keywords: ultrarelativistic 3 He-Au collisions, event-by-event fluctuations, collective flow, femtoscopy 


1. Introduction 

The recent experimental m and theoretical IMS], 
interest in ultrarelativistic heavy-light nuclear collisions 
originates from expectations that studies of such system 
may shed light on mechanisms governing the formation of 
long-range correlations, and thus reveal information on dy¬ 
namics in the earliest phases of the reaction. A successful 
scenario for heavy-light collisions explored in this paper 
involves, exactly as in the well understood case of two 
heavy ion collisions, collective dynamics (hydrodynamics, 
transport models). Collectivity provides in a natural way 
the shape-flow transmutation: the deformation of the ini¬ 
tial configuration (ellipticity, triangularity, etc.) is trans¬ 
formed event-by-event into harmonic flow at freeze-out. 
Moreover, the approximate translational symmetry on the 
initial transverse shape along the spatial rapidity direc¬ 
tion leads to collimated flow, producing the famous ridge 
structures, i.e., azimuthal correlations between particles 
with a large pseudorapidity separation. Yet another vivid 
feature of collectivity is the mass ordering of various ob¬ 
servables mm seen in proton-nucleus collisions m- 

The azimuthal deformation of the fireball in small sys¬ 
tems is due to random fluctuations, as in p-Pb collisions, 
or to a combination of fluctuations and the intrinsic defor¬ 
mation of the small projectile, as in d-Au collisions, where 
the large intrinsic separation between the proton and neu¬ 
tron leads to large elliptic flow, predicted in [5] and ver¬ 
ified experimentally in [3]. Moreover, collisions involving 
projectiles with an intrinsic triangular deformation, such 
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3 He-Au m or 12 C-Au PS], are particularly interesting, as 
they probe systems with nontrivial projectile geometry. 

In this paper we analyze in detail the predictions of 
3+1D viscous hydrodynamics of Ref. m for 3 He-Au colli¬ 
sions at yj snn = 200 GeV, the reaction currently analyzed 
at RHIC. In a previous paper (20] we have proposed spe¬ 
cific tests of collectivity, based on ratios of flow coefficients 
evaluated with 4- and 2-particle cumulants. In the present 
work we focus on other, more direct and immediately ac¬ 
cessible aspects of the reaction, namely, on the elliptic and 
triangular flow coefficients, V 2 and V 3 , and the femtoscopic 
radii. We also explore the ridge formation, both on the 
3 He-side and Au-side in the pseudorapidity direction. A 
comparison to preliminary flow data from the PHENIX 
Collaboration [5T[ indicates a successful description of the 
reaction within our approach. We also confirm the very 
recent findings of Romatschke [22] for the flow coefficients 
and the interferometric radii. 

Throughout this paper we use the three-phase approach, 
consisting of 1) the Glauber [55] Monte Carlo simulations 
of the initial state with GLISSANDO [53] . the interme¬ 
diate 3+ID event-by-event viscous hydrodynamics [TO] , 
and the statistical hadronization at freeze-out simulated 
with THERMINATOR [25]. The three-nucleon config¬ 
urations of the 3 He nucleus are taken from the Green’s 
function Monte Carlo calculations [5B] as provided in [57;. 
The initial entropy density of the fireball in the transverse 
plane follows from the mixed model [281 55i , incorporating 
the wounded nucleons and an admixture of binary col¬ 
lisions, with relative weight a. We use a = 0.125 and 
a NN = 42 mb with a Gaussian wounding profile. The 
entropy density is constructed as a sum of Gaussians in 
the transverse plane of width 0.4 fm, located at the posi- 
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Figure 1: Elliptic and triangular flow coefficients of charged hadrons 
as functions of the transverse momentum, pj-, for 3 He-Au collisions: 
PHENIX preliminary data m (full circles and full squares, respec¬ 
tively) and hydrodynamic calculations (solid and long-dashed lines 
respectively), as well as elliptic flow coefficient for d-Au collisions: 
PHENIX data [321 (empty circles) and hydrodynamic calculations 
(dashed line). 


tions of the wounded nucleons. The assumed longitudinal 
profile of the fireball along the space-time rapidity is dif¬ 
ferent for the left- and right-going participant nucleons; 
we use an approximately linear ansatz for the dependence 
on the spacetime rapidity (see Ref. [3D] for details), that 
fairly well describes the observed asymmetry of the rapid¬ 
ity spectra in d-Au collisions m- 

The centrality class c = 0—5% is selected approximately 
by using events with the number of wounded nucleons 
N VJ > 25. The viscous hydrodynamic evolution is per¬ 
formed with shear viscosity p/s = 0.08 and bulk viscosity 
C/s = 0.04 (for T < 170 MeV). The evolution starts at 
r = 0.6 frn/c, and ends at the freeze-out hypersurface with 
temperature of 150 MeV. The non-equilibrium corrections 
to the local momentum distributions at freeze-out are im¬ 
plemented in the THERMINATOR code. Our analysis is 
carried out with 500 hydrodynamic events, where on top 
of each we generate 5000 THERMINATOR events. 


2. Flow coefficients and long range ridge correla¬ 
tions 

We begin the presentation of our results with the el¬ 
liptic and triangular flow coefficients, evaluated as func¬ 
tions of the transverse momentum, pt , for all charged 
hadrons (Fig. [I]). We evaluate the flow coefficients of 
charged hadrons with \r]\ < 0.5, using the scalar-product 
method [33] El]! with the reference particles taken from 
the Au-side bin, 3.1 < 77 < 3.9 and 0.3 < p± < 5.0 GeV. 
We note large values of both coefficients, in approximate 
agreement with the preliminary results of the PHENIX 
Collaboration ED- For comparison, we also show v 2 for 
d-Au collisions, which is somewhat smaller than in the 
3 He-Au case and compatible with the data. We note that 


Figure 2: Elliptic flow coefficient for charged pions, kaons, and pro- 
tons+antiprotons (solid, long-dashed, and dashed lines respectively) 
from the 3+ID hydrodynamic calculations, plotted as functions of 
the transverse momentum pt- 


a larger value of v<i in 3 He-Au collisions is expected from 
the large initial ellipticity of the fireball [ 20 ] . 

In Fig. [2] we show the mass hierarchy of V 2 (pt), plot¬ 
ting it for charged pions, kaons, and protons+antiprotons. 
The behavior is typical for hydrodynamics, with the p+p 
case significantly lower than pions or kaons at low pt , and 
higher above pt — 1-2 GeV. 

A very vivid manifestation of collectivity is made with 
the help of the two-particle correlation functions, defined 
in the standard way as 


C{Ap, Acp) 
C(A<t>) 


SjA/hMl 

B{At], A</>) ’ 
f dAp S(Ap, Acj)) 
J dAp B(Ap, A(j>) ’ 


( 1 ) 


where the signal S is constructed by histogramming 
the pairs of particles with the relative pseudorapid¬ 
ity Ap and the relative azimuth A (f>, while the back¬ 
ground B is the analogous quantity evaluated with the 
mixed events. The kinematic cuts for the two particles 
correspond to the PHENIX experiment: |?7trig| < 0.35, 
1 GeV < pt ,trig < 3 GeV, and -3.7 < p ass0 c < -3.1 or 
3.1 < i/assoc < 3.9 for the Au-side and 3 He-side, respec¬ 
tively. The result of our model simulation for C(Ap,Acf>) 
are shown in Fig. [3] We note the clear appearance of the 
near- and away-side ridges. 

The result of the projected correlation function C{Ac/>), 
as well as the PHENIX preliminary data [21], are shown 
in Fig. [4] Again, we note clearly the emergence of the 
ridge structure, both in the Au-side and in the 3 He-side, 
in fair agreement with the data. We remark that since the 
kinematic cuts are rather narrow, limiting the number of 
the observed hadrons, it is quite challenging to accumu¬ 
late enough statistics in the hydrodynamic simulation to 
have statistically significant results. Qualitatively similar 
results have been obtained in the AMPT model [31]. 

We note that the precise results for C(Acj)) are sensitive 
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Figure 3: The two-particle correlation function C(A rj. A</>) for a) Au- 
side and b) 'lie-side. The near- and away-side ridges are clearly 
visible. The surface has been smoothed for better visibility. The 
kinematic cuts are specified in the text. 


to the kinematic cuts mimicking the detector acceptance 
and efficiency. Here we have used pt, assoc > 0.1 GeV. 
A higher pr cut would lead to stronger flow, and, conse¬ 
quently, to a higher ridge amplitude. If better accuracy 
is desired, Monte Carlo simulations should carefully incor¬ 
porate the response of the detectors. 

In our THERMINATOR simulations we have incorpo¬ 
rated approximately the transverse momentum conserva¬ 
tion, which is crucial for the reproduction of the features 
of the correlation functions. We have used a procedure 
which accepts only those events where the total transverse 
momentum is sufficiently small. We have verified that 
limiting the acceptance window to \J2iPT,i\ < 5 GeV is 
enough, which corresponds to accepting only about 10% of 
all events {7.. The momentum conservation builds up the 
strength of the away-side ridge, while the near-side ridge 
comes predominantly as a combination of the second and 
third harmonic flow components in the correlation func¬ 
tion. It interesting to note that the observed near-side 
ridge structure for 3 He-side rapidities can be explained 
as an effect of the hydrodynamic expansion of the fireball; 
the collective expansion correlates particles at forward and 
backward rapidities with the event-plane defined by the 
geometry of the fireball. 
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Figure 4: The two-particle correlation function C( A0) for a) Au-side 
and b) 3 He-side. The ridge structure is clearly visible for both cases 
in the data and in the model. The model (empty squares) uses the 
same kinematic cuts as the PHENIX experiment. The preliminary 
data (filled dots) come from Ref. m- 


3. Interferometric radii 

To obtain the Hanbury Brown-Twiss pionic correlation 
radii, we use the two-particle method m to evaluate the 
three dimensional correlation function in the relative pion 
momentum (q Q u t, q s ide, 9iong)- For each given bin of av¬ 
erage transverse momentum of the pair, k±, we construct 
histograms of the correlated and mixed pairs from THER¬ 
MINATOR events [25]. The correlation function, evalu¬ 
ated as a ratio of these histograms, is fitted with Gaus- 
sians, 

C(q,k ± ) = 1+ (2) 

n 2 2 t~) 2 2 t~> 2 2 r> t ~>2 

^o u t4out — A side g side — R long 9 long — 2R out _ long <2out9iong 

The additional mixed term involving q 0 ut®ong appears for 
asymmetric collisions [38] ; including or dropping this term 
does not change the fitted values of the radii R ou t, -Rside> 
and I?iong- Figure [5] presents our results, with stars in¬ 
dicating the 3 He-Au system. We note that the radii fall 
down with k± , as expected from the presence of flow. The 
fall-off is very fast for the small parameter Rout-long- 

The k± dependence of the interferometric radii in central 
d-Au collisions (open circles in Fig. [5]) is similar as for 
3 He-Au, and is in reasonable agreement with the PHENIX 
data [3B]. The 3 He-Au radii are, as expected, somewhat 
larger than for the d-Au case. We note that the ratio of the 
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Figure 5: Pionic interferometric radii -R 0 ut 3 i? s idej ^long? and 
^out—long (Panels a) through d), respectively), plotted as functions 
of the average pion pair momentum for 3 He-Au (stars) and d-Au 
(circles) collisions. The triangles denote the experimental results of 
the PHENIX collaboration for the d-Au system [Hl- 


rms radii of 3 He and d is 1.14, similar to the corresponding 
ratio for the R SI d e radii. 

4. Conclusion 

We reiterate our main results: 

1. The transverse-momentum dependent elliptic and tri¬ 
angular flow coefficients, V2 (j>t) and V3(pt), are as 
in the preliminary data of the PHENIX collabora¬ 
tion [2T] , 

2. The characteristic ridge structure is found, consistent 
with the data [2TJ, both along the Au and the 3 He 
sides. 

3. We have made predictions for the particle-identified 
elliptic flow coefficient V 2 {pt), showing the expected 
mass ordering. 

4. We have evaluated the Hanbury-Brown Twiss inter¬ 
ferometric radii for pion pairs, which are somewhat 
larger than in the d-Au system. They decrease with 
the transverse momentum of the pair. 

All above features are consistent with the collective in¬ 
terpretation of the fireball dynamics formed in ultrarela- 
tivistic nuclear collisions, indicating that evolution at all 
rapidities is driven by the same universal mechanism, de¬ 
termined by the elliptic and triangular deformation of the 
initial fireball, supplemented with fluctuations. 
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